The Cepheid Impostor HD 18391 and its Anonymous Parent 

Cluster 



David G. Turner 1 ' 3 ' 4 , V. V. Kovtyukh 2 , Daniel J. Majaess 1 ' 3,4 
David J. Lane 1 ' 3 , Kathleen E. Moncrieff 1 

turnerOap . smu . ca 



ABSTRACT 

New and existing photometry for the GO la supergiant HD 18391 is analyzed in order to 
confirm the nature of the variablity previously detected in the star, which lies off the hot edge 
of the Cepheid instability strip. Small-amplitude variability at a level of AV = 0.016 ± 0.002 
is indicated, with a period of P = 123 d . 04 ± d .06. A weaker second signal may be present 
at P = 177 d .84 ± d .18 with AV = 0.007 ± 0.002, likely corresponding to fundamental mode 
pulsation if the primary signal represents overtone pulsation (123.04/177.84 = 0.69). The star, 
with a spectroscopic reddening of Es_y = 1.02 ± 0.003, is associated with heavily-reddened 
B-type stars in its immediate vicinity that appear to be outlying members of an anonymous 
young cluster centered ~ 10' to the west and 1661 ± 73 pc distant. The cluster has nuclear and 
coronal radii of r n = 3.5' and R c = 14', respectively, while the parameters for HD 18391 derived 
from membership in the cluster with its outlying B stars are consistent with those implied by 
its Cepheid-like pulsation, provided that it follows the semi-period-luminosity relation expected 
of such objects. Its inferred luminosity as a cluster member is My — —7.76 ± 0.10, its age 
(9 ± 1) x 10 6 years, and its evolutionary mass ~ 19M©. HD 18391 is not a classical Cepheid, yet 
it follows the Cepheid period-luminosity relation closely, much like another Cepheid impostor, 
V810 Cen. 

Subject headings: techniques: photometric — stars: variables: Cepheids — open clusters and associations 



1. Introduction 

The Cepheid instability strip is well recognized 
as the principal region of the H-R diagram occu- 
pied by luminous variable stars of spectral types 
F, G, and K. Less well known is the region blue- 
ward and at higher luminosity than the Cepheid 
instability strip occupied by a variety of stars, typ- 
ically of luminosity classes la to lb, exhibiting vari- 
ability and microvariability at much smaller am- 



plitudes (IMaeder fc Rufenerl Il972l : iRufener et al 
19781 iRufener fe Bartholdilll982l) . Such stars ap- 



pear to obey relationships, as a function of spectral 
type or color, between t heir luminosit i es and semi- 
perio ds of variability (jBurki 1978b; iBurki et al 



19781 ) that, with the inclusion of amplitude as an 



additional parameter, are adequately explained by 
evolutionary m odels of massi ve stars that incorpo- 
rate mass loss (Maeder 198fj|) . 
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The GO la supergiant HD 18391 (jNassau fc Morgan 
has been used previously to 



19521 ) has been used previously to character- 
ize the cool end of the above rela t ionsh ips, tied 
to observations by IRufener et al. I (|l978l) . It is 
listed as a potent ial member of the double clus- 
ter h fc y Pers ei ( Schmidtl 1984 ). as well as Cas 
OB6 (|Stotherdll972h . The star's broad spe ctral 
lines w e re rem arked upon by [Redman (1930) and 
Griffin (1970), although that feature appears to 



1 



be related to high luminosity in late-type stars 
(llmhofJ fl977h . Broad band UBV photo electric 
photo metry for the star was tabulated bv lArguel 
(1966), sased upon 4 observations from Kitt Peak 
National Obse rvatory . Add itional UBV photome- 
try is given by lHauel (|l970h . but the evidence for 



variab ility lies in the observations of Ru fener et al 
( 19781 ). in which the star was observed to in- 



crease gradually in brightness over an interval of 
85 nights during the winter of 1972-73. 

Unlike other stars defining the semi-period- 
lumino sity-color relation for supergiants (jBurkil 
Il978bf ). the semi-period of > 84 days for HD 18391 
was estimated from a single incomplete cycle, and 
is not well established. Our interest in HD 18391 
arose from its spectroscopic parameters, which en- 
abled its interstellar reddening to be estimated us- 
ing a newly-developed technique ( Kovtvuk"rill2007t 
Kovtvukh et al.ll2008h . 

The star's location in the spectroscopic H-R 
diagram (Fig. [J) suggests a possible link to the 
Cepheid instability strip, which raises an intrigu- 
ing question. Does the instability strip for Galac- 
tic Cepheids extend to pulsating GK supergiants 
like those in the Magellanic Clouds, where pul- 
sation periods on the order of 100 days are not 
unusual? The present investigation was initiated 
in order to examine that question in greater detail. 

2. Observations and Air Mass Corrections 

The present study is based upon photomet- 
ric observations of the brightness of HD 18391, 
as well as a detailed s tellar atm osphere study of 



the st ar by iKovtvukhl (|2007l ) and lKovtvukh et al 



(|2008l ) that generated the spectroscopic parame- 
ter T e ff = 5775 K using line ratio techniques, as 
well as logg = 1.2, mctallicity [Fe/H] = 0.02, and 
microturbulent velocity 12.0 km s _1 . New CCD 
F-band differential photometry for HD 18391 was 
obtained on 101 nights between July 2007 and 
April 2009 using the 28-cm Schmidt-Cassegrain 
telescope of the semi-automated Abbey Ridge Ob- 
servatory (ARO) located in Stillw ater Lake, just 



outsi de of Halifax, Nova Scotia (see lMajaess et al 



20081 ). All-sky CCD B V photometry for the same 
field was also obtained on the photometric nights 
of October 22, 23, and 26, 2007, using stars in the 
nearby open cluster NG C 225 as reference stan- 
dards (jHoag et al.l Il96ll h Finally, we obtained 
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Fig. 1.— The location of HD 18391 (numbered filled 
square) relative to other FGK giants and supergiants 
(other filled squares), known Cepheid variables (open 
circles), and the instability strip (dashed lines) in a HR 
diagram that plots surface gravity versus effective tem- 
perature. Evolutionary tracks are depicted for stars of 
4 and 5 Mq, and "BL" denotes the edge of the blue 
loop evolutionary stage. 



CCD spectra in October 2008 for some luminous 
stars in the immediate vicinity of HD 18391, using 
the Cassegrain spectrograph on the Plaskett tele- 
scope of the Dominion Astrophysical Observatory, 
at a dispersion of 60 A mm -1 . 

Differential photometry for HD 18391 was ob- 
tained relative to HD 18327 (AO), with BD+57°668 
(GO: V) serving as a check star. All-sky photom- 
etry for the three stars yielded the results sum- 
marized in Table [21 where comparison estimates 
for the magnitudes and colo rs are taken from mean 
Tycho F-ban d obs ervations dPerrvman et al.lll997t l. 
Argue! (|l966l) . and|Haug] ([l970j). The mean Tycho 
V magnitude for HD 1 8391 wa s adjus ted accord- 
ing to the precepts of iBesselll {2000), but only 
half of the suggested correction to B- V color was 
used, otherwise the result w ould be signif i cantly 
bluer than the values found bv lArgud(|l966l ). lHaugl 
( 197dl ). and this study. The mean Tycho values 
for HD 18327 and BD+57°668, which are much 
bluer, are cited without adjustment. 

Because the comparison stars are much bluer 
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Table 1 

bv magnitudes and colors for program objects 



Star 


V 


B - V 


Source 


HD 18391 


6.937 ± 0.004 


1.966 ± 0.003 


This paper 


(GO la) 


6.936 ± 0.005 


1.982 ± 0.015 


Tycho 1 




6.890 ± . . . 


1.950 ± . . . 


Argue fl966'l 




6.890 ± . . . 


1.930 ± . . . 


Haug (19701 


HD 18327 


8.729 ± 0.010 


0.222 ± 0.013 


This paper 


(AO) 


8.732 ± 0.011 


0.211 ± 0.015 


Tycho 1 


BD+57°668 


9.797 ± 0.001 


0.605 ± 0.007 


This paper 


(GO V:) 


9.827 ± 0.019 


0.676 ± 0.031 


Tycho 1 



Note. — 1 Pcrrvman et al.l U997I) 



than HD 18391, it was necessary to adjust the 
differential photometry for air mass effects, at- 
mospheric extinction being well known to dimin- 
ish the brightness of blue stars by larger amounts 
than that of red stars as air mass increases. HD 
18391 is a yellow sta r of very large reddening 



(jKovtvukh et al. 2008), so its brightness relative 



to bluer comparison stars would appear to increase 
artificially when viewed through large air masses 
if such corrections were not made. Given the low 
level of light variability in HD 18391, tests for 
such an effect are essential to confirm that any 
detected variability originates in the star rather 
than Earth's atmosphere. 

Such an effect can be seen in differential pho- 
tometry bv iFuentes et all (|2006h from Mt. Hop- 
kins for the high latitude hypervelocity blue star 
SDSS J090744.99+024 506.8 (HVS) a t 9 h 07 m 45 s .O, 



+2 o 45'07" (J2000). IFuentes et all (|2006h 



con- 
cluded that the star was a small amplitude vari- 
able because their differential CCD photometry 
relative to redder field stars in its vicinity indi- 
cated a brightness decrease on two nights of ob- 
servation, coincident with an increase in the star's 
air mass. Temporal plots of the brightness of HVS 
on the two nights are replotted in Fig. [2] to dis- 
play its changing brightness as a function of air 
mass X. For comparison stars with colors typi- 
cal of G dwarfs relative to HVS, our experience 
with photometry from Kitt Peak suggests that 
the expected trends should have functional depen- 
dences of about —0.03 for g-band observations and 
—0.01 for r-band observations. The observed de- 
pendences are —0.027 and —0.009, respectively. 



The decreasing brightness of HVS as a function 
of increasing air mass is therefore exactly what 
is predicted for a blue star being compared with 
redder comparison stars, implying that its light 
variations are not intrinsic to the star. 
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Fig. 2. — The dependence of observations by 
IFuentes et all (|2006l ) for SDSS J090744.99 + 024506.8 
on air mass X. Filled symbols represent g-band ob- 
servations, open symbols r-band observations offset 
by —0.1. Straight lines are best fits to the observed 
trends, with residual scatter of ±0.008 in g and ±0.005 
in r. 

The nonvariable status of HVS appears to be 
imp ortant to our u n dersta nding of the star, given 
that Fuentes et al.l ( 20061) based their conclusion 
that it is a slowly-pulsating B-type main-sequence 
star (SPB) primarily on its purported variability. 
Since HVS is not variab le, the origina l spec trum 



of the star obtained by iBrown et al.l (|2005l ) be- 



comes a vital diagnostic tool. The spectrum dis- 
plays broad hydrogen Balmer lines that appear to 
terminate near H13 or H14, which normally cor- 
responds to a su rface gravity o f log<? ~ 4.5 in B- 
type stars (e.g., IVitense 1951 ). That is typical 
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of objects lying below the main sequence, such as 
the blu e horizontal branch stars originally consid- 
ered bv | Brown et al.l |2005t) , and more recently by 
Demarque fc Virata (|2007l ). as a potential origin 
for the object. Perhaps that possibly should be 
reconsidered, given that it is vital to any conclu- 
sions regarding the dynamical origin of hyperve- 
locity stars? 



The Geneva observations for HD 18391 Rufener et al 



(|1978I ). supplemented by a few addi tional unpub - 
lished observations by the group (jBurkil 120071 ). 
were also checked for air mass effects, despite 
their origin from calibrated all-sky photometry. 
The site of the observations is not specified in 
the original references, but was assumed to be 
Gornergrat. The F-band observations are plotted 
in Fig. [3] as a function of inferred air mass, where 
it appears that a small zero-point offset affects 
observations on six of the nights, most of which 
were designated as being of low quality. An ad- 
justment of —0.025 was adopted to bring them 
into agreement with the remaining observations, 
and the resulting data were adjusted for a small 
air mass term in the data (equivalent to an overly 
steep color dependence) and adjusted to a mean 
magnitude near V = 6.94, to be consistent with 
the present data set. 
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Fig. 3. — The dependence of Geneva observations for 
HD 18391 on air mass X. Open symbols represent data 
from six nights of low quality, subsequently offset by 
—0.025 to match observations on other nights. The 
straight line is a best fit to the latter. 

A similar plot is presented in Fig. [5] for our raw 
differential photometry relative to HD 18327. A 
zero-point offset is present midway through our 
data set as a result of a change in CCD detectors, 
which produced a change in the color dependence 
for the observations. There is also a small differ- 
ence in air mass dependence between observations 



in 2007 and 2008-09. Corrections for all such ef- 
fects, an air mass dependence, zero-point offset, 
and differences relative to the results of our all- 
sky photometry, were made to all observations, 
with nightly means summarized in Table [H in- 
cluding a few combined observations from adja- 
cent nights. Plots of the resulting time dependence 
for the Geneva, Tycho, and ARO observations are 
shown in Fig. 
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Fig. 4. — The dependence of Abbey Ridge observa- 
tions for HD 18391 on air mass X. Filled symbols rep- 
resent observations from the first portion of the ob- 
serving period, open symbols observations made after 
a change in CCD detectors. Straight lines are best fits 
to the observed trends. 



3. Period Analysis 

ft is not evident from the available observations 
that HD 18391 is a variable star. The standard 
deviations for the various data sets are ±0.009 for 
Geneva observations, ±0.015 for Tycho observa- 
tions, and ±0.012 for ARO data in Tabled ±0.014 
for the individual observations. The dispersion in 
the ARO check star observations is ±0.016, which 
is larger than the scatter observed for HD 18391 
only because the latter is 3 magnitudes brighter 
than the check star. Its signal therefore had a 
higher signal-to-noise ratio (SNR) than those for 
the comparison stars, and was affected less by pho- 
ton statistics. Small SNRs for the comparison 
stars were inevitable, given that exposure times 
were kept short to avoid saturation of the target. 
The scatter for the check star also appears to be 
random, whereas the observations for HD 18391 
display slight temporal trends. Similar observa- 
tions for other stars observed at the ARO of com- 
parable brightness to HD 18391 exhibit smaller 
dispersion, typically about ±0.008, but the main 
conclusion is that any variability in HD 18391 is 
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Table 2 

F-BAND OBSERVATIONS OF HD 18391 2007-2009 



HJD 




V 


HJD 






V 


HJD 






V 


2454308. 


.737 


6.945 


2454522 


.622 


6 


.935 


2454748 


.599 


6 


.935 


2454309 


.716 


6.943 


2454526 


.589 


6 


.941 


2454750 


.586 


6 


.905 


2454313. 


679 


6.938 


2454532 


.548 


6 


.926 


2454753 


.662 


6 


.924 


2454322. 


605 


6.949 


2454537 


.517 


6 


.920 


2454757 


.658 


6 


.925 


2454324. 


.644 


6.922 


2454539 


.563 


6 


.940 


2454759 


.631 


6 


.936 


2454325. 


638 


6.945 


2454544 


.558 


6 


.933 


2454763 


.616 


6 


.932 


2454332 


.620 


6.944 


2454549 


.589 


6 


.960 


2454764 


.610 


6 


.943 


2454349. 


235 


6.941 


2454551 


.555 


6 


.940 


2454765 


.513 


6 


.946 


2454358 


693 


6.930 


2454559 


.539 


6 


.939 


2454771 


.589 


6 


.941 


2454360. 


776 


6.943 


2454560 


.577 


6 


.942 


2454772 


.643 


6 


.941 


2454362 


661 


6.932 


2454563 


.554 


6 


.934 


2454773 


.613 


6 


.949 


2454364 


.730 


6.925 


2454564 


.526 


6 


.935 


2454783 


.483 


6 


.956 


2454365. 


699 


6.921 


2454565 


.527 


6 


.935 


2454788 


.567 


6 


.950 


2454367. 


.778 


6.910 


2454566 


.532 


6 


.925 


2454804 


.736 


6 


.965 


2454368 


.606 


6.938 


2454571 


.530 


6 


.947 


2454888 


.588 


6 


.941 


2454403. 


461 


6.920 


2454572 


.550 


6 


.940 


2454895 


.579 


6 


.939 


2454412. 


.484 


6.945 


2454573 


.536 


6 


.947 


2454896 


.535 


6 


.947 


2454413 


570 


6.952 


2454578 


.543 


6 


.909 


2454899 


.638 


6 


.929 


2454417. 


.604 


6.954 


2454701 


.568 


6 


.939 


2454901 


.607 


6 


.941 


2454429. 


.452 


6.966 


2454702 


.543 


6 


.948 


2454903 


.570 


6 


.959 


2454433. 


.478 


6.951 


2454714 


.763 


6 


.929 


2454904 


.574 


6 


.958 


2454453. 


496 


6.945 


2454718 


.651 


6 


.937 


2454905 


.636 


6 


.940 


2454456. 


.451 


6.935 


2454725 


732 


6 


.932 


2454906 


.534 


6 


.932 


2454457. 


451 


6.937 


2454727 


.559 


6 


.945 


2454907 


.575 


6 


,944 


2454489. 


.487 


6.934 


2454728 


.744 


6 


.898 


2454908 


.677 


6 


.939 


2454492. 


.680 


6.953 


2454729 


.653 


6 


.929 


2454912 


.574 


6 


,958 


2454497. 


.745 


6.940 


2454730 


,608 


6 


.932 


2454931 


.573 


6 


.949 


2454499. 


.496 


6.937 


2454731 


.667 


6 


.915 


2454936 


.592 


6 


,953 


2454501. 


.545 


6.944 


2454732 


.662 


6 


.938 


2454937 


.577 


6 


.955 


2454516. 


.495 


6.942 


2454734 


.652 


6 


.928 


2454941 


.564 


6 


.952 


2454517. 


.496 


6.944 


2454742 


.777 


6 


.933 


2454946 


.572 


6 


.943 


2454518 


.498 


6.945 


2454746 


.627 


6 


.906 










2454521. 


501 


6.945 


2454747 


.629 


6 


.909 
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at best marginally detectable in existing obser- 
vation s. Previous conclusions by iRufener et al 



(|1978j ) regarding the variability of HD 18391 may 
therefore have been affected by the small air mass 
effect noted earlier. 
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Combin ed ob servat ions fo r HP 18391 
1 19781 ) and iBurkil (|2007f ), top, 



Ruf ener et al . 
iPerrvman et al.l 1 1997 ). middle, and this paper (bot- 



tom) plotted temporally. Open symbols in the bottom 
plot represent observations of the check star for the 
Abbey Ridge observations, offset by —2.35. 

Despite uncertainties about the quality of the 
different data sets, the observations were analyzed 
using various routi nes available in the Peranso 
software package of lVanmunster (l2006h. which in- 
cludes the CLEANest routine (jFoster|ll995|) . Fif- 
teen different period finding routines were used 
with the combined observations, with the average 
value for the most significant period from twelve 
of them being P = 123 d . 04 ± d .06 s.d. A similar 
analysis of observations for the ARO check star, 
BD+57°668, was made as a consistency check, 
yielding a dominant period at P ~ 102 days but 
without producing a realistic light curve. 

The resulting light curve for HD 18391 fits an 
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Fig. 6.— The resulting light curve for HD 18391 
(filled circles) from the adopted ephemeris for P = 
123 d . 04 (upper), along with the best-fitting sine wave. 
Open circles represent the prewhitened observations 
phased to P — 177 d .84 relative to the right-hand scale. 



ephemeris given by: 

HJD max = 2441695.497 + 123.04 E , (1) 

where E is the number of elapsed cycles. It is 
displayed in Fig. [6] The inferred light ampli- 
tude is small, AV = 0.016 ± 0.002, comparable 
to the scatter in the observations and marginally 
significant, and perhaps unusual given that pul- 
sation amplitude typically increases with rate of 
mass loss and luminosity in variable supergiants 
(lMaeder| [l980). The lack of significant light vari- 
ability for HD 18391 presumably implies a small 
rate of mass loss or strong pulsational damping by 
convection. 

A second signal appears in the data at P = 
177 d .84 ± d .18 s.d. after removal of the primary 
signal. The signal is weak, but may be real, given 
that it corresponds with expectations for funda- 
mental mode pulsation if the 123 d signal originates 
from overtone pulsation (123.04/177.84 = 0.69). 
The ephemeris for the second signal is: 



HJD, 



2441735.442 + 177.84 E , (2) 



with a light amplitude of AV" = 0.007±0.002. The 
light curve for the prewhitened observations with 
that period is shown in the lower portion of Fig. [6] 
Conceivably some of the residual scatter in the 
light curves of Fig. [5] results from additional peri- 
odicities in the observations, temporal changes in 
light amplitude, or to rapid period changes in the 
star's pulsation. HD 18391 displays abundance 
patterns typical of stars evolving towards a first 
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crossing of the Cepheid instabil ity strip, includin g: 
a strong lithium feature ( Kovtvukh et al. 2005h . 
As such, radial pulsation in the star should un- 
dergo period_^n^asesof several hours per year 
(e.g., iTurner et al1l2006a ). complicating the anal- 
ysis significantly. The presence of a second pe- 
riod is remarkably similar to what is observed in 
V810 Cen, a F8 la supergiant similar in nature 
to HD 18391 pulsating with periods of 153 and 
104 days a nd displaying tempo r al changes i n light 
amplitude (jKeinzle et al.lll998t lBurkill2006fh The 
light amplitude of HD 18391 is an order of magni- 
tude smaller than that of V810 Cen, so additional 
observations and long-term monitoring are crucial 
to search for similar effects, as well as to confirm 
the small-amplitude variability implied by obser- 
vations to date. 

4. Further Consistency Checks 

Independent checks can be m ade into the na 
ture of the star. For example, iKovtvukh et al 



(2008) derive a spectroscopic reddening for HD 
18391 of E(B-V) = 0.991, but that is based 
upon a different adopted B-V color. With the 
results of Table [21 the reddening becomes E(B- 
V) = 1.020 ± 0.003. A field reddening for HD 
18391 can be found using stars in its immedi- 
ate vicinity. The observed colors o f stars in the 
2MASS catalogue (|Cutri et al.ll2003l ) lying within 
5' of the star are shown in Fig. [71 When ana- 
lyzed using the method illustrated bv lTurner et alj 
(|2008l ). the data imply a field reddening of E(J- 
H) = 0.29 ± 0.03, which corresponds to E(B-V) 
= 1.04 ± 0.09, consistent with the spectroscopic 
result. The obvious A- type main-sequence stars 
in the sample have inferred distances of ~ 650 pc, 
implying that the reddening in the field originates 
relatively nearby, certainly foreground to the GO 
la supergiant. 

The distance to HD 18391 can be estimated 
using luminous stars located in close proxim- 
ity. Table [H lists early-type and luminous stars 
located within ~ 15' of the GO la supergiant 
from Volume I o f Luminous Stars in the North- 
ern Milky Way (jHardorp et al.lll959l ). for which 




photo m etric observ ations are available frq mlHaug 
(jl970l ). Heel (|l970h . and IPerrvman et~aH (|l997h 
Included is the LC class variable M supergiant 
V637 Cas. Spectral types are available for most 



Fig. 7. — A 2MASS color-color diagram, H-K versus 
J-H, for stars within 5' of HD 18391 having magnitude 
uncertainties less than ±0.03. The intrinsic relation 
for main-sequence stars is plotted as a solid line, and 
the same relation reddened by E(J-H) = 0.29 ± 0.03 
as a dashed line surrounded by dotted lines. The kink 
corresponds to AO stars. HD 18391 is depicted by an 
asterisk. 



stars from the ADS data base, including results 
from White fc Wind (|l978l) for V637 Cas, but our 
spectroscopic observations of the B-type stars, 
presented in Fig. [H were used to test and aug- 
ment them. The resulting spectral classifications 
produce the spectroscopic distance moduli given 
in the last column of Table [H where a value of 
R = A V /E{B -V) = 3.2 has been adopted for 
the ratio of total-to- selective abso rption in this 
region of the Galaxy (|Turnerlll97(l . 

All of the luminous stars near HD 18391 are 
reddened by E(B — V) ~ 1.0, with distance mod- 
uli of ~ 11. The exception is the A2 lb su- 
pergiant LS+57°114, which does not fit coevally 
with the remainder of the group. The mean 
unreddened distance modulus for the group is 
Vo — My — 11.12 ± 0.12, corresponding to a 
distance of 1678 ± 96 pc. Interestingly enough, 
there is a sparse group of early-type stars located 
- 10' west of HD 18391 at J2000 co-ordinates 
2:58:28, +57:37:30 (Fig. EJ. The cluster is only 
marginally detectable by eye, although similarly 
sparse groups have been detected previously in 
such fashion by the lead author (e.g., Turnerlfl985t 
Turner et al. 1993). Visually, the cluster displays 
the characteristics of a loose group in the terminal 
stages of dissolution into the field, as confirmed 
by star counts. Ring counts about the adopted 
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Table 3 

Luminous stars within ~ 15' of HD 18391 



Star 


Sp.T. 


V 


HD 18391 


GO la 


6.94 


LS+57°111 


OB+ce,h 


11.42 


LS+57°113 


B0.5 V 


11.03 


LS+57°114 


A2 lb 


11.31 


LS+57°115 


B1.5 Vn 


12.46 


3709 490 1 


B2 Vnn 


11.34 


LS+57°116 


bl v 


11.94 


LS+57°117 


BO V 


10.82 


V637 Cas 


M2 lab-lb 


9.88 



Mean 
s.e. 



- V 


Source 


E(B-V) 


V -M v 


1.97 


1 


1.02 




0.85 


2 






0.81 


1.2 


1.09 


11.13 


1.00 


2 


0.95 


13.18 


0.83 


1.2 


1.08 


11.48 


0.60 


1.4 


0.84 


10.83 


0.89 


2 


1.16 


11.21 


0.87 


1.2 


1.17 


10.77 


2.70 


3 


0.99 


11.31 
11.12 
±0.12 



Note.— Sourc es: 1. This paper, 2. IHauel U9701) . 3. ILeel 1 119701) . 4. 

IPerrvman et al.l 1119971) 
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1.0 


~u 


0.5 - 






M 




ormal 


1.0 


z 


0.5 - 




1.0 




0.5 



TYC 3709-490-1 



400 450 500 

Wavelength (nm) 

Fig. 8. — Low-dispersion spectra for stars near HD 
18391. 



cluster center for J < 15 are shown in Fig. IT0| 
from which one derives nuclear and coronal radii 
of r n — 3.5' (1.7 pc) and R c = 1 4' (6.8 pc), respec- 
tively, in the nomenclature of Kholopov ( 19691 ). 
An unusual feature is the increase in star den- 
sities in the cluster corona, indicative of a ring- 
like feature. A similar characteristic is shared by 
Collinder 70, better known as the Orion Ring, a 
young cluster of stars surrounding e Ori in Orion's 
Belt that is in the advanced stages of dissolution 
into the field. 

The parameters for this group of stars, most 
of which can be identified as early B-type dwarfs 



from their 2MASS colors (Fig. fTTj) . are similar 
to those for stars lying closer to HD 18391. The 
implied reddening is E(J-H) — 0.40 ± 0.02, corre- 
sponding to E(B-V) = 1.43 ± 0.07, and the dis- 
tance modulus is Jo — Mj = 11.07 ± 0.15, cor- 
responding to a distance of 1638 ± 113 pc. The 
adopted mean of the two estimates is d — 1661±73 
pc (11.10±0.10 in intrinsic distance modulus). HD 
18391 and most of the luminous stars of Table 2] 
lie in the cluster corona, a by-now commo n char- 
acteris tic of many young star clusters (e.g., Burkil 
1978a]). The predicted number of cluster mem- 



bers for J < 15 lying within 3.5' of the cluster 
center is 16 ± 8 according to the star counts, com- 
pared with exactly 16 identified in this region from 
2MASS photometry. The total number of cluster 
members with J < 15 according to the star counts 
is 102 ± 35, typical of a poorly-populated cluster. 

If HD 18391 is a member of the cluster and 
B-star group, its implied luminosity is My = 
—7.76 ± 0.10. It appears to fit in an evolution- 
ary sense as well, according to the unreddened 
color-magnitude diagram of Fig. [12] for stars in 
the B-star group. The stars appear to have an 
age of log* = 6.95 ± 0.05, or (9 ± 1) x 10 6 
ye ars, according to an isochrone fit to models 
bv iMevnet et al.l (|l993l ). The irregular M super- 
giant variable V637 Cas does not match the same 
isochrone very well, almost certainly because of 
heavy mass loss during it s previous evolution (e.g., 
Jura fc Kleinmann 19901) . A similar situation ap- 



plies to the M supergiant BC Cyg a s a member of 
Berkeley 87 ([Turner fc Forbes! 1 19821 : iTurner et af 
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Table 4 
Radial Velocity Data. 



Star 


JD 


V R (km s- 1 ) 


HD 18391 


2449739.3755 


-41.3 


HD 18391 


2452711.3171 


-44.4 


HD 18391 


2453785.2410 


-43.5 




Fig. 9. — The view of the anonymous cluster located 
at 2:58:28, +57:37:30 (J2000) on the Palomar Obser- 
vatory Sky Survey blue plate. The field measures 30' 
by 30' with north up, and HD 18391 is the bright star 
on the left edge of the field. [The National Geographic 
Society-Palomar Observatory Sky Atlas (POSS-1) was 
made by the California Institute of Technology with 
grants from the National Geographic Society.] 



2006ri iMaiaess et "all |2008h. and t he M3 lab-lb 
member of NGC 2439 (|Turne ij[l977h . The mass of 
stars at the red turnoff (RTO) for a 9 x 10 6 year- 



old cluster is 18 \M { 



© 



according to iMevnet et al 



(|1993l ). implying that the mass of HD 18391 is 
~ 19M©. Possible membership of HD 18391 in 
Cas OB6 can be ruled out. Cas OB6 members 



inclu de evolved 04 and 05 stars (jHillwig et al 



2006), implying a much young er age, and are more 



distant: 2.01 kpc according to Mcl'nik & Efremov 



1995), ~ 1.9 kpc as adopted by iHillwig et al 



20061 ). They also have radial velocities of ~ —50 
km s _1 , marginally more negative than the mea- 
sured velocities for HD 18391 summarized in Ta- 
ble [4] Attempts were made to measure radial ve- 




15' 20' 25' 30' 

Ring Radius 

Fig. 10. — Ring counts for J < 15, with their Poisson 
uncertainties, for the anonymous cluster centered at 
2:58:28, +57:37:30 (J2000). The field star level based 
upon rings with r > 13' is shown, as is a possible 
trend for cluster stars. The dimensional parameters 
are r n = 3.5' and R c = 14'. 



locities for the associated B stars observed spec- 
troscopically in Fig. but with negative results 
because of zero-point problems. The three coolest 
stars in the sample were found to have nearly iden- 
tical velocities, however. 

The implied luminosity of HD 18391, M v = 
— 7.76, or Afboi = —7.79, agrees closely with expec- 
tations from published semi-period- luminosity re - 
lations for Cepheid-like supergiants (|Burkilll978bl h 
provided that the longer period (P = 178 d ) ap- 
plies. Interestingly enough, both HD 18391 and 
V810 Cen also fit predictions from the Cepheid 
period-luminosity relation, as shown in Fig. 1131 
The results for the two FG class la supergiants are 
based upon their assumed me mbership in A non 

lTurnerlll98j) . al- 



(HD 18391) and Stock 14 (see 
though the case for V810 Cen as a cluster member 
has been questioned on the basis of a conflict with 



the r adial velocity of Stock 14 stars (jKeinzle et al 
Il998l ). The case for HD 18391 also requires radial 
velocity confirmation. 

The radius of HD 18391 implied by its luminos- 
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-0.5r 

h-k ; 

0.0 - 



HD 18391 Cluster 




10 - 

11 - 

12 - 

13 - 
14- 
15 - 




1.0 1.5 
J-H 



Fig. 11. — A 2MASS color-color diagram, H-K ver- 
sus J-H (upper section), for stars within 5' of 2:58:28, 
+57:37:30 (J2000) having magnitude uncertainties less 
than ±0.05. The intrinsic relation for main-sequence 
stars is plotted as a solid line, as well as for E(J-H) 
= 0.40. Filled circles represent likely cluster members, 
open circles probable field stars. Note the large num- 
ber of reddened stars blueward of the AO star kink. 
The lower section is a color-magnitude diagram for the 
same stars, indicating the best zero-age main-sequence 
fit for J-Mj = 11.07 ± 0.15. 



ity and effective temperature is 329i?0, whereas a 
value of 583i? Q would ch aracterize a Cepheid o f 
the same pulsation period ( Turner &; Burkel l2002) . 
The implied radius also produces a value of log g = 
0.69 with the effective temperature inferred from 



the m odel atmosphere analysis ([Kovtvukh et al 



2008), consistent with the cited value of log g =1.2 
when non-LTE and possible mass-loss effects are 
taken into account. But HD 18391 is clearly not 
a Cepheid. Its variability presumably originates 
from pulsation, radial pulsation if the implied pe- 
riods are correct, highly damped by convection. 
Variability possibly linked to the rotation of HD 
18391 would imply a rotation rate of ~ 90 km s _1 , 
much too large for a star that has expanded by 
two orders of magnitude since the main sequence 
stage. The color variations of the star are also 




1.0 1.5 
(B"V) 

Fig. 12. — An extinction-free H-R diagram for neigh- 
boring stars to HD 18391, the circled point. Isochrones 
are shown for logt = 6.9 (upper) and \ogt = 7.0 
(lower) . 



1978bD . further constraining t he nature of pulsa- 
tion in the star (Maeder 1 19801) . 




1.5 2.0 
log P 



2.5 



in antiphase to the brightness variations (jBurki 



Fig. 13. — The period-luminosity relation for Galac- 
tic Cepheids associated with open clusters and associa- 
tions. Classical Cepheids are denoted by filled circles, 
and the Cepheid-like supergiants V810 Cen and HD 
18391 by open circles. 



5. Discussion 

The suspected, long-term, low-amplitude, light 
variability of HD 18391 is confirmed by recent 
observations of the star, although the light curve 
derived from recent CCD observations contains 
considerable scatter, a problem that may be as- 
sociated with the difficulty of making accurate 
extinction corrections for observations from a 
site located near sea level. The dominant pe- 
riod of variability of P = 123 d . 04 is reasonably 
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solid, but implies a semi-period of 61 d .52, smaller 
than the value of > 84 d originally estimated by 



Rufener et al.l (|1978l ). The latter value lies closer 
to P = 177 d .84 (semi-period of 88 d .92) found 
for the weaker second signal. A luminosity of 
My = -7.76 is implied for HD 18391 from likely 
membership in a surrounding group of luminous 
stars that, with the GO la supergiant, are outliers 
to an anonymous, poorly-populated open cluster 
~ 10' to the west that appears to be dissolv- 
ing into the field. The likely age of the B-star 
subgroup surrounding the GO la supergiant is 
log* = 6.95 ± 0.05, or (9 ± 1) x 10 6 years. The 
luminosity estimate is consistent w ith the existin 
semi-period-luminosity relation of iBurkil (1978 
for such objects. 

It is interesting to note that both V810 Cen and 
HD 18391 fit the period-luminosity relation for 
cluster and association Cepheids extremely well, 
in fact unerringly in the case of HD 18391, de- 
spite the fact that their implied effective temper- 
atures and luminosities place them well to the hot 
side of the classical Cepheid instability strip. The 
mechanism responsible for the variability of such 
Cepheid impostors is still poorly constrained, al- 
though non-radial pulsation and effects tied to 
mass loss have been considered previously ([Maeder 
19801 ). The present results for HD 18391 suggest 
the possibility of radial pulsation. It is curious 
how well Cepheid-like supergiants obey the classi- 
cal period-luminosity relation for Cepheids. It is 
not clear why that should be, although their pres- 
ence among the luminous population of galaxies 
may be important to the identification and use of 
Cepheids as distance indicators, particularly for 
low-amplitude variables. 
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